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and is contrasted with that of the Aig phonon. A mechanism involving the charge 
transfer fluctuation between the two oxygen ions in the Cu02 plane coupled to the 
crystal field perpendicular to the plane is discussed and the resulting electron-phonon 
coupling is evaluated. Depending on the symmetry of the phonon the weight of 
different parts of the Fermi surface in the coupling is different. This provides the 
opportunity to obtain information on the superconducting gap function at certain 
parts of the Fermi surface. The lineshape of the phonon is then analyzed in detail 
both in the normal and superconducting states. The Fano lineshape is calculated 
in the normal state and the change of the linewidth with temperature below Tc 
is investigated for a d^2_y2 pairing symmetry. Excellent agreement is obtained for 
the Big phonon lineshape in YBa2Cu307. These experiments, however, can not 
distinguish between ci^.2_y2 and a highly anisotropic s-wave pairing. 
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I. INTRODUCTION 



As a result of the recent dispute as to whether the superconducting gap in the high-Tc 
materials is s— or (i— type all experiments which can contribute to resolve that question 
are of special importance It was shown recently that the study of the polarization 
dependence of the electronic Raman scattering can address these questions and can distin- 
guish between the different orientations {dx2-y2 or dxy) of the gap as well. This evidence 
can be made more complete by including the phonons and the effect of the electron-phonon 
coupling, which is the subject of the present paper 

The analysis of the lineshape of optical phonons in the high Tc superconductors has 
been the focus of a large body of investigation Of the many Raman active modes in 

tetragonal superconductors, all the modes that have been observed in the cuprate materials 
using in-plane polarizations transform according the Aig symmetry except one, which trans- 
forms according to Big symmetry [Q. Since this phonon obeys its own selection rules, the 
frequency of the mode can be unambiguously assigned, which is in contrast with the early 
confusion of the assignments of the separate Aig modes. 

In this paper the phonons involving the perpendicular vibrations of the in-plane oxygens 
are investigated. The displacements of the two oxygens 0(2) and 0(3) are in phase for Aig 
symmetry and out of phase for Big as is shown on Fig. (1) for the case of YBa2Cu30y. The 
coupling of these phonons to the electrons in the plane has an interesting feature, namely a 
completely flat Cu02 plane taken out of the surrounding atomic environment has the mirror 
symmetry through the plane. Thus the linear coupling between these phonons and the in- 
plane electrons is absent. In a crystal however, as it has first been pointed out by Barisic 
et.al. 0, that symmetry can be broken since there is an electric field perpendicular to the 
plane due to the surrounding ions forming an asymmetric environment. This is in contrast 
to the La compound, where only the small tilting of the octahedra breaks the symmetry ^ . 
The perpendicular electric field in the 1:2:3 material can be responsible for the buckling of 
the Cu02 plane, i.e. the slight separation of the planes formed by the Cu and O ions. This 
small distortion however cannot be responsible for the large electron-phonon coupling to be 
discussed here. 

Since the long-range charge transfer fluctuations between unit cells are screened, only 
intracell charge transfers must be considered. Let us denote the charge transfer fluctuations 
at 0(2) and 0(3) by 6px and 6py (see Fig. 1.) respectively. A charge transfer fluctuation 
where 6p-j: — 6py 7^ is obviously coupled to the electric field and the crystal deforms in 
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the form of the Big phonon, similarly Spx + Spy 7^ is coupled to the Aig phonon. The 
fluctuations of the in-plane charges are either due to the transfer (i) between the in-plane 
oxygens and other ions, or (ii) between the two in-plane oxygens. The fluctuations of the 
Aig symmetry 6px = Spy can be a consequence of the transfers between the in-plane copper 
and oxygens. For Big symmetry both (i) and (ii) may be realized. 

Considering the first possibility the charge transfer in 1:2:3 between the bridging 
oxygen 0(4) and the 0(3) is a good candidate. As long as the non-bonding orbital of 0(4) 
parallel to the CuO chain has a partial electronic density of states near the Fermi surface, 
the transfer breaks the x — y symmetry. That mechanism, which is sensitive to the position 
of the electronic partial density of states of the bridging oxygen relative to the Fermi energy, 
has been worked out in detail [Q, but the recently observed Raman spectra are not consistent 
with the predicted d^y superconducting gap 0. 

The second possibility does not depend on such details of the electronic band structure, 
since only the conduction band plays a role. The charge transfer fluctuation 6px + Spy = 
between 0(2) and 0(3) involves especially those parts of the cylindrical Fermi surface which 
are close to the and ky axis. Thus the influence of the electron-phonon coupling on 
the electronic Raman continuum depends strongly on whether a large superconducting gap 
opens at that part of the Fermi surface. More precisely, the electronic contribution can 
be characterized by the azimuthal quantum number m on the cylindrical Fermi surface 
0J^, and m = ±2 corresponds to Big symmetry, while out of the two (m = 0, m = ±4) 
channels of Aig symmetry only m = ±4 is relevant since the m = mode describes charge 
fluctuations between cells, which are screened by the long-range Coulomb interaction. The 
different Fermi surface areas probed by Raman scattering in the Aig and Big geometry are 
illustrated in Fig. (2). {Aig and Big transforms like cos(40) and cos(20) respectively). 
Thus assuming a d^^^yi type superconducting gap the Big mode probes the areas with the 
largest gap, while the Aig averages the areas of both the largest gap and the nodes. The 
most appropriate tool in the Raman experiments to probe the electron-phonon coupling is 
to study the Fano interference which is formed due to the simultaneous scattering of light 
on the phonons and the electronic continuum (see eg. Refs. @] and [§). The fit of the Fano 
lineshape provides information not only on the coupling but also on the electronic continuum 
influenced by the opening of the superconducting gap. The appearance of a sharper gap in 
the Big symmetry than in the Aig provides a unique identification of the gap of d^-i^y-^ type 
0. It must be emphasized however, that these experiments are not sensitive to the sign of 
the order parameter, therefore can not distinguish between ci-wave superconductivity and a 
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very highly anisotropic s type. 

The present status of the relation between experiment and theory can be summarized 
as follows. The vibration of the Big mode which appears at roughly 340 wavenumbers in 
all cuprate materials, is connected with the antisymmetric out-of-plane vibrations of the 
0(2) and 0(3) ions in the Cu-0 plane 0. The net charge transfer of this vibration in 
the unit cell is zero and thus the long range Coulomb forces are incapable of screening 
the charge fluctuations. Consequently, the mode has a large cross section for scattering 
incoming photons and thus appears in Raman experiments as a large sharp signal centered 
at frequency shifts corresponding to the energy of the mode. Due to the strength and 
unambiguous identification of the mode, the Big phonon has been lavished with attention 0. 
In particular, the spectral lineshape of the mode has been a subject of intense scrutiny both 
in the normal and superconducting state of the cuprates. In the normal state, the asymmetry 
of the lineshape indicates strong mixing of the phonon with the electronic continuum and 
fits to a Breit-Wigner or Fano lineshape have been made with great success This has 
usually led to the assertion that the electron-phonon coupling at least for this mode is large. 
However no estimates of the coupling from a microscopic theory have been presented except 
for the theoretical suggestion by Barisic and Kupcic 0. Further, the anomalous changes 
of the Big phonon as the material enters into the superconducting state have also been 
exhaustively documented P,|5|,p!0[| . The temperature dependence of the changes in the Big 



phonon's lineshape have been used to help determine the magnitude of the energy gap of the 
superconductor [|1^]. Recently, a theory of Raman scattering in (i— wave superconductors has 
been presented and while the general features of the theory fit well with the experiment, 
no detailed fit of the lineshape could be made without further knowledge of the mechanism 
and strength of the electron-phonon coupling. 

In the present paper we provide a detailed theory for the behavior of the Big phonon in 
both the normal and superconducting states of the cuprate materials using the mechanism 
suggested by Barisic et.al. 0J^. In particular, we investigate the mechanism of electron- 
phonon coupling resulting from crystal field effects and describe the Fano lineshape in the 
normal state. The theory is then generalized to the superconducting state and the lineshape 
is calculated in detail for a superconductor with d^2_y2 pairing. 

The paper is organized as follows. In Sec. II. we develope the mechanism which leads 
to first order electron-phonon coupling due to the presence of the crystal field. Based on 
the 3-band model for the Cu02 plane the electron-phonon coupling constant is evaluated. 
In Sec. III. we apply these results in order to fit the experimental data in the normal state. 
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Particular attention is paid to the Fano resonance. Section IV. is devoted to the behavior 
of the phonon hneshape in the superconducting state. Here we show that the temperature 
dependence of the spectrum is due to the change in the electronic response, and the data 
are consistent with a dx2-y2 gap. Finally our conclusions are given in Sec. V. 



II. MECHANISM 



In the first part of this section we define our model for the electrons in the CUO2 plane. 
This allows us to introduce notations used in the second part, where we develop a microscopic 
theory of the coupling of these electrons to the out-of-plane phonons due to the crystal field 
perpendicular to the plane. 



A. The model 



Using the notations of Ref. |Tl| we consider the 3-band model for the CUO2 plane de- 



scribed by the following Hamiltonian: 

^° = ^EOn,<x + ^ E + (1) 

n,0,(T 

where bl^ ^ creates an electron with spin a at a copper lattice site n, while ^ annihilates 
an electron at one of the neighboring oxygen sites n + 6/2 determined by the unit vector 
6 assuming the four values, (±1,0) and (0,±1). An oxygen atom between the two copper 
atoms at sites n and n + (5 is labeled by either (n, 5) or (n + 6, —6). Moreover, e = Ed — Ep 
is the difference of the Cu and O site energies, t is the Cu-0 hopping integral and = ±1 
depending on whether the orbitals (with real wavefunctions) have the same or opposite sign 
at the overlap region. Assuming Cu dx2-y2 and O p orbitals = —P§, and we can choose 
-P(i,o) = 1 and -P(o,i) = —1- For simplicity, we neglect direct 0-0 hopping. 
The momentum representation is defined by the following formulae: 

Ee^''""^k,., (2) 



and 



,5,,, = ^ E exp[ika(n + 5/2)]a„,k,^, (3) 
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where a is the lattice constant, N is the number of Cu sites and a is x or y ior S — (±1, 0) 
and S — (0, ±1) respectively. In this representation the Hamiltonian in Eq.(l) decouples for 
different momenta as 



(4) 



k,<T 



where 



Hk,a = + {i6^[a^,k,atx(k) - ay,k,atj/(k)] + h.c.}, 



and 



taiia) = 2tsm{aka/2). 
H^^ is then easily diagonalized as 

/3 



(5) 



(6) 



(7) 



where /? assumes the values +, — , and for the antibonding, bonding, and nonbonding 
bands, respectively. The corresponding electronic energies are ii^o(k) = and 



with 



£;±(k) = e/2 ± ^(£/2)2 + Q2(k)^ 



(8) 



(9) 



The transformation between the original (a, b) and the new (d) electronic operators is de- 
scribed by 



( ^ \ 

ax 



(e+,eo,e- 



do 



(10) 



where the column vectors of the transformation matrix are given by 



ty 



(11) 



and 



e± 



( \ 

-itx 



(12) 



In the last three equations we have dropped the k, a indices for clarity. 

In the physically relevant situation for the Cu02 plane the upper band is close to half 
filled. Since the three bands do not overlap in the present model, we will consider only 
a reduced one-band Hamiltonian describing the upper band and neglect the bonding and 
nonbonding bands. Then the reduced Hamiltonian is given by 



(13) 



k,(T 



with 



E{k) = e/2 + ^(£/2)2 + 02(k). 



(14) 



Furthermore, the expressions of the transformation from d operators to a and b operators 

reduce to 



(15a) 



where 



(15b) 
(15c) 



,(k)=E(k)[E2(k)+l]2(k)]-V2^ 



(16a) 



0,(k) = -it,{\i)[E\\i) + Q^(k)]-V2, (16b) 

0^(k) = ity{y)[E\\i) + l]2(k)]-V2. (16c) 

As we will see in the next subsection, this simple model contains the most important ingre- 
dients in order to describe the electron-phonon coupling. 
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B. Electron-phonon coupling 



As we have discussed already, the electrons in the Cu02 plane do not couple in first order 
to the phonon modes with displacement vectors perpendicular to the plane, because in that 
case the hopping integrals change only in second order in the displacements. However, if 
there is an electric field perpendicular to the plane, first order electron-phonon coupling is 
generated. In our case this field originates from the asymmetric environment around the 
Cu02 plane (see Fig.l.), and is also responsible for the buckhng of the plane with restoring 
force provided by the covalent bonds. 

Let us consider the effect of a lattice periodic crystal field E(r) = — V^extlr) on the 
electrons in the plane. The electron density at each (displaced) site couples to the external 
field via the Hamiltonian 

H' = -eEn,<7{^n,aV<T0ext[«n + Uf,(an)] 
+ «n,x,cT«n,x,<70ext [on + ax/2 + (an)] 

+ 4,y,a«n,y,a</>ext[an + ay/2 + Uj/(an)]}, (17) 

where u;,(an), u^(an), and Uj^(an) arc displacement vectors of the Cu, 0(2), and 0(3) in 
the unit cell at the lattice site n, e is the electron charge, and x and y are unit vectors in 
the corresponding directions. Expansion in the displacements up to first order leads to 

H' = H'^ite + Hel-ph + (18) 

where H'^^g renormalizes the copper-oxygen site energy difference e only, while the term 
linear in u generates an electron-phonon interaction 

+ E3;Ua;(an)a^_^an,x,o- 

+ EyUy{aii)ai^y^„an,y,a}, (19) 

where E^, Ej,, and Fiy are the electric fields at the Cu, 0(2), and 0(3) sites respectively, 
which are parallel to the z-axis due to symmetry. In the absence of chains — 'Ey would 
hold. 

Hei-ph can be written in momentum representation with the help of Eqs.(2-3), and 
Eq.(15) allows us to express the interaction of phonons and electrons of the reduced Hamil- 
tonian Eq.(13) in the usual form 
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H, 



el—ph 



^ EE 9\iK q)4+q,A,<x[cA(q) + 4(-q)]- 



(20) 



q,A k,(T 



Here C\{q) annihilates a phonon mode A with wavevector q and gx(k,c[) is the couphng 
constant of that mode to an electron of wavevector k. Based on our microscopic model the 
coupling constant can easily be evaluated using standard procedures of the quantum theory 



of phonons |]12 



In the context of Raman scattering we are interested in q = optical phonons. In case 
of the Big phonon the Cu displacement is zero, while the 0(2) and 0(3) atoms have equal 
and opposite displacements (Fig.l.b.). The corresponding coupling constant is evaluated as 



£/Bi„(k,q = 0) = 



n 



\ 2MoUJB,, V2 



(21) 



where Mq is the oxygen mass and oje^^ is the phonon frequency. In the following we make 
the approximation {Eix)z = (Ey)z = E^. Due to the opposite displacements of the oxygens, 
the difference of the 0q, functions (given by Eq.(16)) determines the k dependence of the Big 
coupling. Since Q(k) = const, on the Fermi surface (see Eq.(14)), this coupling constant is 
proportional to t^.(k) — t^(k) oc cos(aA;^) — cos(aA;j^), which is clearly of Big character. In fact, 
this is exactly the second order Fermi surface harmonic of our model band structure 
(Eq.(14)) transforming according to the Big symmetry. Therefore we can write the coupling 
in the form 



where 



\cos(akr) — cosfaA;,, 



(22) 



(23) 



is the normalized second order Fermi surface harmonic of Big symmetry with the normal- 
ization constant at the Fermi energy Ep given by [T^ 



Nb^Ep) 



J cPk6[EF - i?(k)][cos(aA;^) - cos(aA;y)]^ 



(24) 



J d^k5[EF - E{k)] 

All information about the strength of the coupling is now compressed in the expansion 
coefficient 

' 2^ I Nb,,{Ep) 



9b 



Is 



h 



\ 2MoUJB,, Ef{2Ef - e) 



(25) 
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In case of the Aig phonon the oxygen displacements are the same both in magnitude and 
in direction (Fig.l.a) and the displacement of the copper is negligible due to the relative 
rigidity of its vertical bonds. The same procedure employed for the Big phonon yields the 
following coupling constant for the Aig phonon 



^Ai„(k,q= 0) = eE, 



^ ' [\M^)f + (26) 



In this case however, the coupling depends on k only through the energy -E'(k), i.e. it is 
constant on the Fermi surface for any band filling. Therefore within the present model the 
Aig phonon couples to homogeneous density fluctuations only, and since these fluctuations 
are suppressed by the long range Coulomb interaction, the resulting coupling is vanishingly 
small. This conclusion remains valid if we allow for a finite Cu displacement as well |]14| , 
but does not necessarily hold if e.g. 0-0 hopping is included. 

At the end of this section we wish to evaluate the coupling in Eq.(25). For a half 
filled band the Fermi energy Ep = e/2 + ^ i^/^y + (2t)^, and the normalization constant 
NBj^giEp) = 4. It is reasonable to suppose that e/2 <^ 2t, and in this limit we only need 
the electric field Ez at the oxygen site for a numerical value of gsig- In order to estimate 
this field we suppose that the charges on the planes of the unit cell of YBa2Cu307 are 
evenly distributed, and the charges (per unit cell, in units of e) for the planes Y, Ba, and 
Cu02 are +3, +2, and -2 respectively, while the remaining (chain) region has -3 electron 
charge to insure neutrality of the unit cell. Each plane produces an electric field at the 
oxygen site we are concerned with independent of its distance. Since all planes in the 
sample contribute, we consider ever larger environments of the Cu02 plane in question. We 
calculate the field produced by the two neighboring (Y and Ba) planes first, then include 
the next nearest neighbors etc. The series of values for E^ is of course not convergent, but 
has a period of the unit cell. Although a more accurate calculation applying the Ewald 
summation is desirable 0, for an estimate we use the average value in this series, which 
yields eE^ = —2'ne^lc? = — 6.1eV/A. According to Eq.(25) this crystal field generates a 
coupling gsj^g = 0.12eV. The relevant dimensionless coupling A^i^ = ^rg's^Jft'^Big = 0.078 
if we use a total density of states at the Fermi level Np = 0.22/eV corresponding to a 
bandwidth of 9eV. We note here that due to the presence of chains the tetragonal symmetry 
of the Cu02 plane is weakly broken, and (Ex)z ^ (Ey)^. This leads to the slight mixing of 
the Aig and Big modes. 
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III. NORMAL STATE 



In this section we utilize our calculations for the coupling constant to discuss the resulting 
spectra of the cuprates in the normal state, while in the next section we discuss the changes 
in the spectra due to superconductivity. 

The inelastic scattering of photons from a metal can either be caused by collisions with 
phonons or via the creation of electron-hole pairs. As well, the phonons interact with 
the electronic continuum and have a dynamical effect on the way photons are scattered. 
The total light scattering cross section resulting from these contributions is depicted by 
the Feynman diagrams in Fig. 3. The coupling constant 7 describes the electron-photon 
coupling, (i.e., the Raman vertex for particle-hole creation due to the vector potential of the 
incoming light). In the limit of small momentum transfer and frequency transfers smaller 
than the optical band gap, the Raman vertex is simply related to the curvature of the band 
dispersion e(k) and the incident (e^) and scattered (e^) polarization light vectors via 

which can be expanded in terms of Fermi surface harmonics 0™'(k), 7(k) = J2mlm4>^0^)- 
The remaining vertices we denote by Qp-p for the photon-phonon vertex, and g\ for the 
electron-phonon vertex as discussed in the previous section. The bare phonon propagator is 
given by 

where is the intrinsic phonon linewidth resulting from e.g. the decay of the phonon caused 
by the presence of an anharmonic lattice potential. In the metallic state, the electron-phonon 
coupling renormalizes the phonon propagator such that 

D^{uj) = ^^l^""] , ^ , (29) 



where x\ is the complex electronic susceptibility evaluated in channel A [T^. Defining the 
renormalized phonon frequency through 

Qjl = ^l{\-\{yj)), (30) 

where = 2g1x'x{^)/^\ and x' denotes the real part of the susceptibility, we arrive at 
the following expression for the renormalized phonon propagator 
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where Vx is the total frequency dependent hnewidth of the phonon, = V\ + glxx{uj). 

We now sum up the diagrams in Fig. (3) for the Raman response in terms of the 
susceptibihty xx- After some lengthy but trivial algebra we obtain 

\2 ( 



II ( \ _ [^ + ^a) 12 /// N 



+%%rU [1 + A(^)//?]4, (32) 

\LjJ -\- UJa, 



with 



/? = ^^^, 7A=(0A(k)7(k)). (33) 

Here (■ ■ ■) denotes an average over the Fermi surface. Thus 7a represents the symmetry ele- 
ments of the Raman vertex projected out by the incoming and outgoing photon polarization 
vectors. 

Eq. (32) is a generalized form of the Breit-Wigner or Fano lineshape describing the inter- 
action of a discrete excitation (phonon) with an electronic continuum. Here the frequency 
UJa = cjaVI + 1^ sets the position of the antiresonance of the lineshape. If the intrinsic 
phonon width r\ = then at the antiresonance Xxjuiii^ = ^a) = 0. Since xa describes 
the electronic contribution to Raman scattering, which is generally featureless and a smooth 
function of frequency, we can replace everywhere xx by the value it takes on near the phonon 
frequency to fit the Raman spectra in the vicinity of the phonon. However to fit the entire 
spectrum in the normal state, the full uj dependence is required. 

The description of the continuum in high Tc superconductors has been the focus of a 
large body of both theoretical and experimental work. While at present no theory can 
completely describe the full symmetry- dependent continuum for a large range of frequency 
shifts, we now give expressions for the continuum in two separate cases. In the presence of 



impurity scattering, the continuum lineshape has a simple Lorenztian form ||T6 



UJT 



X';.{uj) = Np—-^. (34) 



UJ^ + Tx 



Here = l/rm=o — ^/tx is the channel-dependent impurity scattering rate reduced by 
vertex corrections. While this form can fit the data on the cuprate materials at least in the 
Big channel (i.e., crossed polarization orientations aligned 45 degrees with respect to the 
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a — b axis in the Cu-0 planes) at low frequencies, the high frequency tail which remains 
relatively constant and the frequency behavior at low frequencies in other channels cannot 
be accounted for. 

Another form for the continuum can be arrived at by taking nesting of the Fermi surface 
into account as has been done by one of the present authors . The Raman susceptibility 



has a similar form as Eq. (34) with the impurity scattering rate replaced by the addition of 
the electron-electron and impurity scattering rates on a nested Fermi surface and the effect 
of mass renormalization taken into account, 



oj turn* (tu) /m, 

2a 

with m*{uj)/m = 1 H In 

TT 



y'(/3'T)2 + ^2 



(35) 



where and uOf. are constants determined by a fit to the electronic continuum in the 
normal state. This form for the Raman susceptibility provides an adequate description to 
the continuum in the normal state of both YBCO and BSCCO (see Ref. [0) and thus we 
will use this form rather than Eq. (34). 

In order to fit the Raman spectra in the normal state, one sees from Eqs. (32) and 
(35) that a great deal of parameters are required. A fitting procedure is now discussed 
that greatly reduces the number of free parameters. The first step in the procedure is to 
determine the bare phonon parameters ujx and F^, which in principle can simply be read off 
from a fit to the lineshape of the phonon in the insulating state where g\ = While the 
intrinsic width of the phonon can be read off as F^^ = 2.5cm~^, a problem arises with the 
position of the Big phonon, since it shows only a small renormalization with doping from 



the insulator to the metallic state |^ , the origin of which is unknown. This could be due in 
part to changes of the lattice parameters with doping and / or the tetragonal to orthorhombic 
transition |^. Therefore, we thus must keep this parameter free and use the bare frequency 
we obtain from fitting the normal state. While in principle Qp-p can be determined as well 
at this stage, since the intensity of the spectra is given in arbitrary units, this parameter 
(which sets the overall intensity) must also remain free. 

Turning next to the metallic state, the continuum parameters a, and 7a are tuned to 
fit the full frequency range of the continuum minus the contribution of the phonons. We then 
tune the renormalized phonon position, c^Big, which we later can compare with the result 
predicted from the value of gsig obtained in Section II. Thus the only unknown parameters 
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are Qp-p (and thus which also sets the antiresonance position tUa, and the bare phonon 
frequency. 

With these constrainsts, such a procedure produces the theory hues shown in Fig. (4a) 
compared to the data of Ref. |]19| on YBCO. The parameters used to obtain the fit are as 
follows: ujb^, = 358cm-\ ub,, = 348 cm-\ Tb.^ = 4.9 cm-\ T^^^ = 2.5 cm-\ a = 0.55, (3' = 
3.3, qI^^Nf/ujb,, = 0.0624, r^7^^ = 600 cm~\ujc = 12, 000cm-\and^p_p/7B,^ = 0.0026, 
which sets Ua = 360 cm~^. Similar fits have been obtained before via the usual Fano 
expression p0|j5[] albeit with unconstrained parameters determined solely via a fitting routine. 
The value of the coupling constant used is roughly 20 percent smaller than the one given by 
the crude approximation in Section II. Given the uncertainty involved in the band parameters 
and the neglect of screening (which will reduce the crystal field and thus lower the coupling 
constant), we remark that the coupling constant derived from a microscopic picture of 
charge transfer within the unit cell provides an accurate description of the strength of the 
asymmetric Fano lineshape of the Big phonon in YBCO. 

IV. SUPERCONDUCTING STATE 

Turning now to the superconducting state, we note that Eq. (32) is still a valid de- 
scription of the Raman lineshape provided that we use the Raman susceptibility calculated 
for the superconducting case. In fact, it is possible to obtain important information about 
the symmetry of the electron pairing in a superconductor through an investigation of the 
changes in the phonon lineshape below Tc. In particular, the electron-phonon coupling de- 
pends crucially on the symmetries of both the phonon and the order parameter. The role 
of the coupling between a phonon of given symmetry and an order parameter of a different 
symmetry was recently discussed by one of the authors , where it was shown that changes 
in the phonon lineshape are the greatest for a phonon which possesses the same symme- 
try as the order parameter while smaller changes are predicted for phonons of a symmetry 
orthogonal to that of the energy gap. This is pictorially shown in Fig. 2, which shows 
the energy gap A(k) squared and the electron-phonon vertex (yf(k) squared, which enters 
into the Raman response function x as a weighted average around the Fermi surface. Thus 
the phonon vertex and the energy gap will contructively (destructively) interfere with each 
other if they have the same (orthogonal) symmetry. This fact can be used to determine the 
predominant energy gap symmetry. 

In Ref. [0] the electronic contribution to Raman scattering xx was calculated in detail for 
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a superconductor with dx2^y2 pairing symmetry and good fits were obtained to the Raman 
spectra of BSCCO. For details of the theory, we direct the reader to that Reference and 
simply write down the Raman response function obtained for the Big channel using a d^^-y-i 
gap, A(k, T) = Ao(T) cos(2(y9) for a cylindrical Fermi surface (here x = cu/Aq): 

Sttx \x[(l + 2x2)ir(l/x) -2(l + x2)E(l/x)], X>1. ^ ^ 

The real parts are obtainable through Kramers-Kronig analysis (see Ref. j^]). The resulting 
response has a logarithmic divergence at the pair edge (an artifact of the two-dimensional 
nature of the Fermi surface) and vanishes as for small frequencies. To make a fit to 
the continuum in the superconducting state, we first convolute Eq. (36) with a Gaussian 
to mimic the effect of finite z dispersion of the Fermi surface, experimental resolution, etc. 
Fitting the resulting function to the superconducting state continuum for YBCO leads to 
value A = 240 cm~^ and a smearing width F/A = 0.2 (not to be confused with the intrinsic 
phonon linewidth). Then, Eq. (36) can be used to draw a fit to the full response using 
the parameters obtained from the fit to the response in the normal metal. 

The resulting fit to the lineshape using Eq. (36) and its corresponding real part (deter- 
mined via Kramer's Kronig transformation) at T = 20i^ is shown in Fig. (4b) compared to 
the data on YBCO from Ref. ||2l]]. Here the parameters used are the same as in the normal 
metal case (Section III) with the only change resulting in the use of Eq. (36) and its real 
part rather than Eqs. (34-35). Using the coupling constant A = 0.0624 as obtained in fitting 
the normal state leads to cj^ig = 337 cm~^ and F^^^ = 7.9 cm~^, i.e., the phonon broadens 
and softens at low temperatures compared to its normal state lineshape. 

The fit points towards the predominance of an energy gap of Big symmetry for the 
following reasons. First we remark that the theory predicts the correct magnitude and sign 
of the phonon renormalization due to superconductivity. This would be in marked contrast 
if the gap was s— wave or of d— wave nature but of other symmetry {d^y, d3z2_j.2,etc). In the 
first case an isotropic gap would lead to trivial coupling of the energy gap to the symmetry 
of the phonon and all phonons would show the same qualitative behavior, which is not the 
case for the cuprates 0. In the second case, previous work has demonstrated that the 
channel probed via Raman scattering which shows a peak at the highest frequency in the 
continuum in the superconducting state of all the channels gives the predominant symmetry 
of the energy gap. This again points to d^^^yi pairing. Further, if the gap was of another 
symmetry (i.e., or Aig) then the position of the phonon, cj^^^/A = 1.42, would predict 
small phonon softening since (i) the real part of the phonon self energy changes from negative 
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to positive values at roughly uj/A = 1.5, and (ii) the weight (7^(k) | A(k) p) over the Fermi 
surface which determines the coupling would be smaller for a gap of different symmetry than 
the vertex. Therefore, while although Raman cannot distinguish between and energy gap 
which changes sign or not around the Fermi surface, the evidence from the data indicate 
that the electronic pairing is predominantly of d^^-y^ symmetry. 

V. CONCLUSIONS 

In order to improve the study of the superconducting gap anisotropy in high temperature 
superconductors the interactions of the electronic continuum and the characteristic 
phonons shown in Fig. 1 have been included in the theory. These interactions and the 
intrinsic width of the phonons had complicated a previous comparison of the theory for the 
continuum for certain gap anisotropy to the experimental data since the contribution of the 
phonons had to be subtracted in an intuitive way (see Ref. 2). In addition, the phonon 
lineshape is of course distorted by the Fano interference. 

The starting point of our calculation was the model calculation of the electron-phonon 
interaction based on Ref. 6 and a simple tight binding calculation for the electrons is 
used. The latter is certainly not very adequate but demonstrates the correct characteristic 
symmetry properties of the electronic wave functions. 

Our main conclusions are as follows: 

(i) Phonon lineshape - The asymmetric Fano lineshape of the Big phonon is reproduced 
while the interaction of the continuum with the Aig phonon is much weaker according to 
the experiments and disappears in our model calculation (see Sec. II B). Using more 
adequate electronic wave functions, this coupling may reappear but still must be weaker 
than in the case of the Big phonon; 

(ii) Electric Field - The perpendicular electric field acting on the oxygen atoms in the 
Cu02 plane has been determined by a comparison with the experiments in the normal 
phase, assuming tetragonal symmetry, i.e., E = = Ey. This symmetry is certainly 
destroyed by the presence of the chains as shown in Fig. 1, and results in a mixing of the 
Aig and Big symmetry channels. It is obvious from the coupling presented in Sec. II B 
and from the calculation of the diagrams depicted in Fig. 3, that starting with, e.g., light 
polarizations assigned to Big symmetry, the coupling with E^ 7^ Ey leads to a contribution 
to the electronic continumm of Aig symmetry which is further coupled to the Aig phonons. 
The Aig phonons do appear in the Big spectrum observed by experiments (see Fig. 4a), 
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but it is not clear whether this is dominantly due to misahgnments of the samples or to the 
channel mixing. A more elaborate calculation of the electric fields and Ey arising from 
different ions would certainly help to resolve this question; 

(iii) Gap anisotropy - The excellent fit of the combined electronic continuum + phonon 
spectrum in Fig. 4b is obtained by assuming the simplest form of the superconducting gap of 
da;2_y2 type for Big symmetry. The parameters are borrowed from the fitting of the normal 
state. The only additional parameters introduced are the gap amplitude defined above Eq. 
(36) and an additional smearing F ~ 0.2Ao with a non-unique origin (see Sec. IV). This fit 
certainly supports the previous fit in Ref. |^ where the phonons had been subtracted. 

Concerning the anisotropy of the gap it has already been pointed out that Raman scat- 
tering provides information only concerning the absolute value of the gap. The fit of the 
low energy spectra in the present paper clearly demonstrates that assuming d— wave pairing 
in the superconductor, the gap with dx2-y2 symmetry gives an excellent fit but d^y, d^z, and 
dyz are not applicable (see e.g., Ref. [0]). A very highly anisotropic extended s— wave state 
which has predominantly Big symmetry cannot be ruled out. 

For instance recently the photoemission experiment of Ref. suggests a gap which can 
be represented by the form A{ip) = Aq[1 + a cos(4y9)]/(l + a), with the parameter a just 
exceeding one. Our fits clearly show that a fit with a = 1 are not satisfactory. The reason 
is that the different contributions of the different gap features to the low energy (tu << Aq) 
part of the spectra with e.g. Big symmetry: (i)the large gap around ip = ^n, (ra = 0, 1, 2, ■ ■ ■) 
does not contribute (see Fig. 2); (ii) the gap changing sign at = ^ + contributes as 
uj^ at low frequencies; (iii) the gap touching zero with zero slope (parameter a = 1) at 
ip = ^ + contributes as u at low frequencies (see Ref. [Q]); (iv) the gap with zero slope 
but not at zero energy (parameter a ^ 1) gives a sharp discontinuity and approximately 
linear term (modulo logs) above or below it (a < 1 and a > 1, respectively). 

On the basis of the above feautures, the zero slope at energy larger than ~ 20cm~^ and 
at (f = \ + can be ruled out. Any feature with a gap minimum on a smaller energy 
scale can however still be possible, but the region of the Fermi surface where the minima 
occur must be very narrow compared with those given by the expression mentioned above. 
As the minima should occur in an anomalously narrow region, the slopes coming or leaving 
the minima must be very large. Such a behavior appear to us as very unlikely. 

In summary, the main result of the present paper is the confirmation of the conclusion 
of Refs. , namely, that the Raman spectra are in excellent agreement with a supercon- 
ducting gap of dj.2_y2 symmetry, but gaps with very sharp features on a small part of the 
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Fermi surface cannot be ruled out. 
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FIGURES 

FIG. 1. The unit cell of the Cu02 plane in YBa2Cu30y is shown with the atomic displacements 
corresponding to the phonons of Aig (a.) and Big (b.) symmetry. A finite electric field E 
perpendicular to the planes due to the asymmetric environment (Y above, Ba below) is indicated. 
This induces the respective charge transfer fluctuations {6p) denoted by the two way arrows. 

FIG. 2. The dashed line above the cylindrical Fermi surface (F.S.) represents a dx-^_y2 super- 
conducting gap |A(i?i))p ~ cos^(20). The areas shaded by the vertical dotted lines mark the regions 
of the F.S. which are probed by Raman scattering in the Aig (a.) and in the Big (b.) geometries. 

FIG. 3. Feynman diagrams depicting the three contributions to Raman scattering in metals. 
Here the solid line is the electron Green's function, the wavy line the phonon propagator, and the 
dashed line the photon propagator. The vertices are defined as in the text. 

FIG. 4. (a). Fit of Eq. (32), (34-35) to the Big spectra from Reference ll on YBCO at 
T = IQQK. The parameters used to obtain the fit are discussed in the text;(b) Fit of Eqs. (32) 
and (36) to the Big spectra from Reference on YBCO at T = 20i^. The parameters used to 
obtain the fit are discussed in the text. The additional phonons besides the 340 cm~^ have been 
subtracted off for clarity (see Fig. 4a). 
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